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Pion dissociation by the Mott effect in quark plasma is described within the generalized Beth- 
Uhlenbeck approach on the basis of the PNJL model which allows for a unified description of bound, 
resonant and scattering states. As a first approximation, we utilize the Breit-Wigner ansatz for the 
spectral function and clarify its relation to the complex mass pole solution of the pion Bethe-Salpeter 
equation. Application of the Levinson theorem proves that describing the pion Mott dissociation 
"^^q" ■ solely by means of spectral broadening of the pion bound state beyond Tmom leaves out a significant 

aspect. Thus we acknowledge the importance of the continuum of scattering states and show its 
, role for the thermodynamics of pion dissociation. 
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I. INTRODUCTION 



We investigate the thermodynamics of mesonic bound states in hot medium consisting of a nonidcal quark plasma 
with correlations in the pion and sigma meson interaction channels. A special emphasis is put on a correct description 
of the bound states' dissociation in the vicinity and beyond the Mott temperature. To this end we utilize the 
Polyakov-loop-extended Nambu-Jona-Lasinio model at finite temperature. 

The model has two order parameters: the chiral condensate, determining the value of the dynamically generated 
quark mass m(T) related to the chiral symmetry breaking/restoration transition, and the Polyakov loop variable 
<£>(T) which is an order parameter for deconfinement. Their values are obtained from a self-consistent solutions of the 
coupled gap equations for these parameters which correspond to the location of the minimum of the thermodynamical 
potential Q(m, $;T) in the mean-field approximation (MFA). Going beyond the MFA, we evaluate the contributions 
*~Q , from pion and sigma meson fluctuations within the Gaussian approximation to the path-integral representation of the 
■ thermodynamic potential. Above the Mott temperature one observes the spectral broadening of bound states leading 
ly-j | to the appearance of complex mass poles in the mesonic propagators. A proper analysis leads to coupled Bethe-Salpeter 
equations from which meson masses and corresponding spectral widths are obtained. Next, the equation of state for 
the quark-meson system is obtained in form of a generalized Beth-Uhlenbeck equation [lj which describes the effects 
^v^j . that chiral symmetry restoration and deconfinement have on the contributions from pions and sigma mesons and 
y—{ ' allows for an adequate description of bound, resonant and scattering states on equal footing. The spectral functions 
in the mesonic channels are evaluated from the complex-valued polarization loop integrals. An examination of the 
Levinson theorem proves that in order to correctly describe the Mott effect, the continnum of mesonic correlations 
(scattering states) has to be taken into the account. As a result of this analysis, we obtain a description of pion 
dissociation where the effect of the vanishing bound state is exactly compensated by the occurrence of a resonance in 
the continuum of scattering states. 

The present work improves on previous works within the NJL model [H, where unphysical quark degrees of 
freedom appeared in the hadronic phase due to the lack of the coupling to the Polyakov-loop and the Stefan-Boltzmann 
limit was not obtained due to a misplaced momentum cutoff. We also improve the recent work [B| by properly discussing 
the role of Levinson's theorem in the PNJL model approach to the generalized Beth-Uhlenbeck equation of state for 
the quark-meson plasma. The model presented joins both exact limits of finite-temperature QCD - the pion gas at low 
temperatures and the quark-gluon Stefan-Boltzmann limit at high temperatures - within a microscopic chiral quark 
model approach. 
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II. QUARKS AND LIGHT MESONS IN THE PNJL MODEL 



The definition of the Nambu-Jona-Lasinio model [1]-[12[ with Polyakov loop [l3|j-[17] is given in the Appendix [Al 
Here we start from the expansion of the thermodynamic potential 



il(T t ri = -~hiZ[T,V,n] 



(1) 



around homogenous mean field values, which leads to the decomposition of the auxiliary fields a' and ir' into their 
mean field and fluctuation parts 



cr' = cr M F + c, 7?' = 7? (ttmf = 0) 



(2) 



Utilizing the decomposed fields results in the factorization of the partition function into the respective parts describing 
mean field and fluctuation contributions: 



Zmf[T, V,fx] = cxp 
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Here the mean field inverse propagator is 

s mf = 7° fan - M + A ) -j-p-rriQ - a MF = 7°(iw„ - fi + A ) - j-p-m 



(4) 



(•5) 



The thermodynamic potential in the mean field approximation of the PNJL model is evaluated to be given by the 
following expression [l7| 
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where factors Nf, N c originate from performing the trace operation and are a consequence of isospin and color 



symmetry. The quark energy is given by E p — \/p 2 + m 2 , are defined as E^ = E p =f fx, and 
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(7) 
(8) 



In the mean field approximation of the PNJL model the values of the constituent quark mass m and the Polyakov- 
loop variable $, along with its complex conjugate $, are obtained from the condition that the thermodynamic potential 
should be minimized with respect to these parameters, which is augmented by the stability conditions. For /j, = 0, we 
have <3? = $ and thus the minimizing conditions are given by 
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MF 
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These conditions are equivalent to a set of coupled gap equations flU . For the mass gap equation we get 



m = m + AN f N c Gs 



d 3 p m 
Wf 3 E~ P 



l-U(E p )-f+(E p ) 



where 



f*(E p ) = [$ e -^ E ^ +2$e-W E ^ +e- 30 ( E ^]/N*(E p ) 



(9) 



(10) 



(11) 



are the so-called generalized Fermi functions, characteristic for the PNJL model. One should note that if $ — > 1, the 
expression pip reduces to the standard NJL model Fermi functions. For PNJL calculations we should additionally 
find the values of $ from corresponding gap equation [17| at given T and fi. 
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FIG. 1: Temperature dependence of the quark masses m(T) and Polyakov loop variable 4>(T) at /i = 0. Results for quark 
mass are scaled by m(0) = 367.5 MeV. 



TABLE I: The set of model parameters reproducing observable quantities (in brackets) and {qq) 1 ^ 3 = —240.772 MeV [l8l ]: 



m [MeV] 


A [MeV] 


G S A 2 


F-K [MeV] 


M-x [MeV] 


5.495 


602.472 


2.318 


(92.4) 


(135) 



In order to solve (|10[) . a set of model parameters has to be determined: the cutoff parameter A, the current quark 
mass too (in the chiral limit too = 0) and the coupling constant Gs- These parameters arc fitted at T = to reproduce 
physical quantities: the pion mass M n = 135 MeV, the pion decay constant F v = 92.4 MeV and the quark condensate 
(qq) 1 /3 = -240.772 MeV. The used parameters [lj] are shown in Table |U 

Since NJL-type models are non-renormalizable, it is necessary to introduce a regularization, e.g., by a cutoff A 
in the momentum integration. Following [17], in this study, in case of thermodynamic quantities, we use the three- 
dimensional momentum cutoff for vacuum terms and extend the integration till infinity for finite temperatures. A 
comprehensive study of the differences between the two regularization procedures (with and without cutoff on the 
quark momentum states at finite temperature) has been performed in [19]. 

Solutions of the gap equation (fTU)) and the corresponding gap equation for the Polyakov loop variable $ at zero 
chemical potential and nonzero T are presented in Fig. [T] Above the critical temperature, which is equal T c = 237 
MeV in the chiral limit and T c = 251 MeV for a finite current quark mass, one observes chiral symmetry restoration 
indicated by the rapid decrease of the constituent quark mass and the Polyakov loop variable becoming close to 1. The 
corresponding mean field contribution to the pressure, given by P = —£Imf> is shown on Fig. |TJ] (scaled by a factor 
T 4 ) in section VI. One observes that in opposition to classical NJL models (e.g., @, 5), quark degrees of freedom arc 
suppresed below the critical temperature in the PNJL model. Moreover, above the critical temperature, the gluonic 
degrees of freedom are accounted for correctly. 

The contribution to the thermodynamics stemming from the fluctuations described by Q is evaluated in a scheme 
where we expand the logarithm up to the second (Gaussian) order according to 

°° x k 1 

\n(l-x) = -Y^ T = -x~-x 2 + ..., \x\<l, (12) 
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to get 



7(2) 

J FL 
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(13) 



where we have introduced E[<7, if] = a + 175x7?. 

Performing the calculation leads to the subsequent factorization of the thermodynamic potential into parts describ- 
ing the contribution from mesonic correlations corresponding to a and if channels of interaction 
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with the polarization loop IIm {qo, q) given explicitly by 
n M (q ,q) = -N C N^ 
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One easily obtains the thermodynamic potential corresponding to the chosen meson part of the partition function up 
to the Gaussian order given by 







(16) 



where the degeneracy factor cLm equals 1 for sigma mesons and 2 for pions. 

From the point of view of the polarization operators, the pseudoscalar and scalar meson masses can be defined by 
the condition that for q 2 = M\ , th e corresponding polarization operator YIm(M^) leads to a bound state pole in 
the meson correlation function pjj]. For mesons at rest (q = 0) these conditions correspond to the Bethe-Salpeter 
equations 
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Solutions of the two Bethe-Salpeter equations (flT]) and (fl~8|) constitute the set of meson masses and are presented 
on Fig. [5] (left panel). The Mott temperature, which is defined by the condition M n (Ty[ 0ti ) = 2m q (Ty[ ott ) , is for 
given parameters TMott — 231 MeV in the chiral limit and TMott — 256 MeV away from it. The modification of 
quasiparticle properties is clearly visible: up to the Mott temperature TMott, the a mass practically follows the 
behaviour of 2m q {T), with a drop towards the pion mass, signalling chiral symmetry restoration (in the chiral limit 
the a mass exactly coincides with twice the quark mass up to TMott)- In the same region, the pion mass remains 
practically constant (and equal zero in the chiral limit). At T ~ TMott however, the masses of chiral partners become 
approximately degenerate, M a ss M n , and then both masses increase linearly with temperature. 

The corresponding meson pressure is shown on Fig. [2] (right panel). 



III. GENERALIZED BETH-UHLENBECK APPROACH TO THE QUARK-MESON PLASMA 

Deriving the thermodynamic equation of state for a quark-meson system in the Beth and Uhlenbeck form reduces 
to the introduction of scattering phase shifts into the formula for the thermodynamic potential. This requires that 
we analytically continue the propagator into the complex plane. 

In a first step, we introduce meson spectral functions A 9 M (uj,q) by utilizing the integral representation of the 
logarithm in (| 16[) . followed by the usual expression of the propagator by means of A 9 M (u>, q) [20| . i.e., 

= - d9^-i —, T = - / + °° I " &A 9 M (u,$. (19) 

Jo 2 5 2 _L_n M ( eo ,9) J-00 2nq -LjJ 2g* mV ' q > V 1 

Using the fact that the spectral density is given by the discontinuity of the propagator at the real axis we arrive at 

[ GS d 9 A9 , * ■ [ Gs d 9 f Q g , , , Q9 , ■ *\ { l-2G s n M (u-iri,q) \ 

L ® = ~ l J V [ Sm + Ql SM q) ) = ~ l { l-2GsU M (u J + i r ] ,q ) ) ' (20) 
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FIG. 2: Temperature dependence of the meson masses Mm (left panel) and meson pressure for 7r (solid lines) and a (dotted 
lines). 



where the argument of the logarithm is by definition the scattering matrix Sm(w, q) in the Jost representation [3j. This 
normalized complex function can also be represented by means of a scattering phase shift S>m(^, q) — exp[2i$^/(o;, q}], 
which allows us to identify 



f Gs da 

/ -^A 3 M (w,q) = -HnS M (oJ,q) = 2<f> M (oj,q) . 
Jo z 9 

Utilizing the above identity leads, after performing the Trace operation, to the following result 
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(22) 



where Dm(s) is a generalized mass distribution (density of states) containing all the dynamics of the system. The 
above expression is the so-called generalized Beth-Uhlenbeck form of the thermodynamic potential. The connection 
to the original Beth-Uhlenbeck expression for the second virial coefficient can be found in [3[ in detail. 

For an analysis of meson masses and corresponding widths, necessary to obtain a description of mesonic correlations 
as bound states dissolving into resonant states above the Mott temperature, one splits the polarization function 
IIm^Oj 0) according to 

(23) 



(24) 



n M [qo, 0J = 4N c N f h - 2N c N f P M I 2 {q ) = h - P M h(qo), 
where, in the limit q = that we utilize, the integrals ii and I 2 are given by 



h = 
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2E P — qo 



2E P + q 



(25) 



with Pm = —qo for pions and Pm = —qo + 4m 2 for sigma mesons. By identifying qo = Mm — i\F m one can perform 
the complex mass pole analysis leading to 
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AN c NfG s 
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(Re I 2 {M M ) - 1 ImI 2 (M M )) , 



(26) 
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which decomposes into coupled Bethe-Salpeter equations for meson mass and meson spectral width. 

In the first departure beyond the pole approximation (Appendix [Bjl , it is justified above TMott to consider D^j(s) 
to be described by a Breit-Wigner type function 



A R (s,T) = a R 



M M T 



M 



M 2 M ) 



(M M T M ) 



(27) 



where Mm is the meson pole mass, V m is the corresponding meson width and a R is a normalization factor. Below 
TMott, where the spectral broadening T(T) of the states vanishes, the above expression becomes the delta function 
typical for the spectral function of a mesonic bound state. The meson phase shift $m corresponding to (|27[) should 
be of the form 



$m(s) » (j) R {s) 



arctan 



( 4m 2 -Ml, \ 
\ M M r M J 



arctan 



s M M 



M M T 



M-L M 



— arctan 



Am 2 - M M 



M M T 



M 1 M 



(28) 



where by introducing the notation 4>r(s) we acknowledge the fact that the above phase shift is connected with resonant 
properties of the mesonic correlations. 



IV. LEVINSON'S THEOREM FOR QUARK-MESON THERMODYNAMICS 



In order to inspect the validity of the approach so far presented, in our analysis we consider the Levinson theorem 

d$M 



ds 

4m 2 ^ 



U7T , 



(29) 



where n denotes the number of bound states below the threshold 4m 2 . Indeed, it is easy to check that the resonant 
phase shift <pji(s) alone does not fulfill (121)1) . This implies that the scattering phase shift should be composed of at 
least two parts. In fact, as was demonstrated in it is appropriate to decompose the scattering phase shift 'I'm into 
a part corresponding to the mesonic correlation and a part describing quark- antiquark scattering, 



$M =<t>R + 4>sc ■ 

Namely, using (j2"Tj) we can represent the total scattering phase shift $m as 



M 



— In ■ 



2 l-2G s IL M (uj-ir],q) ' 
Then it is straightforward, using (|23[) and the relation between logarithm and arctan functions, to show that 

2G s P M lmI 2 



$m = — arctan 



and, by several more consecutive steps, to obtain 



1 - 2G s h + 2G s P M Re h 



(30) 



(31) 
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At this point it is enough to recognize the above expression as the sum of arctans to finally obtain 
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which proves the formula pop , where 
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^Re7 2 / 



(35) 
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and 



t)R = arctan 



1G s h 



Im/o 



2G S \I 2 



Pi 



M 



l-2Gsh 
2Gs\l2\ 2 



Reh 



(36) 



Using the conditions (|26l) . in the above equation masses and widths could be identified in accordance with 

As the first part of this decomposition is independent of mesonic properties, we presume that it is connected with 
the scattering states' input to the thermodynamics. On the other hand, the second part of (f3"4"|) describes solely the 
behavior of dissolving mesons. Thus we are assured that the correct description of mesonic correlations accounts not 
only for bound and resonant states' contribution, but also for the input from the scattering states. 

In our analysis we will use a combined approach, where the scattering part of the phase shift <&m is defined according 
to (1341 and the resonant part is given by a delta function below Tmou and by the Breit-Wigner ansatz (l27l) beyond 
it. Explicitly, we take 



1 # M (s) 

D M {s) = j 

7T as 



(5(s 



Ml) 



T < T r 



Mott 



(37) 



T > Tmou 



We will in the following regard the states for positive real s > 4m 2 only. 

Finally, the scattering states' contribution to the density of states is normalized and exactly compensates the 
contribution from the resonance (or the bound state, resp.). This accordance with the Levinson theorem confirms the 
validity of the presented approach. 



V. RESUTS AND DISCUSSION 



The above described complex mass pole analysis leads to the results for the meson masses Mm and widths Tm as 
presented on Fig. [3l We observe that the sigma meson width is non-zero for all the temperatures considered, although 
below Tiviott it is not significant and therefore allows to consider sigma to be a quasi bound state. Above the Mott 
temperature pion and sigma masses quickly become equal and so do their spectral widths. This is the imprint of the 
chiral symmetry restoration where the a and 7r meson, being chiral partners, become degenerate. 
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In what follows we will concentrate on the discussion of the tt meson, because it undergoes the Mott transition from 
the bound state to the resonant correlation in the continuum, accompanied by a jump of the scattering phase shift at 
threshold from 7r to in accordance with the Levinson theorem. 

In Fig. 2] we show the phase shift $m (lower panel) in the decomposition (I30|) into its resonant (upper panel) and 
scattering continuum (middle panel) parts, obtained from the solution of equations (pj|) and (pJ5)) . 

The pion pressure resulting from utilizing the Breit-Wigner spectral function alone and together with the input 
from the scattering states is shown on the Fig. [5] along with the pressure of the massive pion gas. Fig. [6] shows 
the pion pressure as compared with the quark pressure obtained from the mean field approximation of the PNJL 
thermodynamic potential. We notice that including the description of the scattering states results in a correct 
description of pion thermodynamics in the vicinity and beyond the Mott temperature. 



VI. CONCLUSIONS 



Pion dissocation by the Mott effect in quark plasma was described within the generalized Beth-Uhlenbeck approach 
on the basis of a PNJL model whereby a unified description of bound, resonant and scattering states was given. 
As a first approximation, we utilized the Breit-Wigner ansatz for the spectral function and clarified its relation to 
the complex mass-pole solution of the pion Bethe-Salpeter equation. It has been demonstrated that a description of 
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FIG. 6: Comparison of the pion pressure obtained with the input from the MFA approximation to thermodynamic potential. 



the pion Mott dissociation solely by a spectral broadening of the pion bound state when it enters the continuum of 
unbound states for temperatures beyond TMott necessarily entails a violation of Levinson's theorem. In order to solve 
this problem we have extended the approach beyond the complex mass pole approximation and solved the scattering 
phases in the pion channel of quark-antiquark interaction. The account for the scattering continuum in accordance 
with the Levinson theorem leads to a strong reduction of the pion pressure above the Mott dissociation temperature. 
We suggest that the behavior of the scatterinng phase shift in the pion channel and its temperature dependence across 
the Mott transition, as obtained in the present work, can be used to develop a generic ansatz for the behavior of 
hadronic densities of states to be used in a generalized PNJL-hadron resonance gas model [2l| that embodies the 
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Mott dissociation of hadrons. 
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Appendix A: Nambu— Jona-Lasinio model with Polyakov loop 



Confinement in pure SU(N C ) gauge theory can be simulated by introducing an effective potential for a complex 
Polyakov loop held. The PNJL Lagrangian [l3j]-[l7| is 



Cpnjl = q (*7^ - m - 7 V) Q + ° M ^ T Mqf ~ U($[A],&[A];T) . 



(Al) 



The quark fields are coupled to the gauge field A^ through the covariant derivative = 9„ — iA^ . The gauge field 
is A^ — 6q A — i6%A4 (the Polyakov gauge). The field $ is determined by the trace of the Polyakov loop L(x) [H| 



= — Tr c L(x) 



(A2) 



where L(x) = "Pexp 



dr Ai (x, r) 



Tju are the vertices for the scalar (a 1 ) and pseudoscalar (7?') four- 



fermion interaction channels. The gauge sector of the Lagrangian density (|A1[) is described by an effective potential 
U ($L4],¥L4];T) fitted to the lattice QCD simulation results in pure SU{3) gauge theory at finite T [TH [la] with 



U (<I>, $; T) _ b 2 (T)- b. 



rp4 



b 2 (T) = a + ax ( — 



$$ - -1 4- <I>" j + — (<[><!> 
To 



—3\ &4 /— tN 2 
1 



T, 



, 



T, 



, 

"■■> l y 



(A3) 
(A4) 
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TABLE II: Parameters of the effective potential WL4]. 



Oo 


Ol 


a 2 


as 


b 3 


bi 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 



The parameters of the effective potential (|A3I) and (|A4[) are summarized in Table HU 

In general, the parameter T depends on the number of active flavors and the chemical potential. In the present 
work we use T = 208 MeV as has been proposed in [22[ . 

The partition function in the path-integral representation is then given by 



Z PN j L [T,V,fi] = ypTfPgcxpj^ dT J v d3x 
+ G s (qT*>qf + G s (qT„, q y - U($[A],$[A];T) 



(A5) 



where the interaction vertices are written explicitly. By means of the Hubbard-Stratonovich transformation we are 
able to ingrate out the quark degrees of freedom to arrive at the partition function written solely in terms of collective 
fields 





'_ 

















4G S 



U{$[A];T) 



Tr In 



where S 1 [a' , if'] is the inverse propagator given explicitly as 



S x [a' = 7°(iw„ - n + A ) - 7 • p- m - a'Y a > - n'T w > . 
and the operation Tr is taken over color, flavor, Dirac and momentum indices of quark fields. 



(A7) 



Appendix B: The Breit-Wigner ansatz 

We introduce the following pion propagator 

S M (P 2 ) 



1 



l-2GTI(p 2 ) ' 

In the first step we expand the polarization function around the mass pole 



no? 2 ) = n(p 2 = A/ 2 ) + {p 2 - M 2 ) 



and we insert this expansion into the defined propagator 



dn{ P 2 ) 

dp 2 



+ . ■ 



(p2=M 2 ) 



Sm(p 



1 



all 



'HI 



1 - 2GII0 2 = M 2 ) -2G(p 2 - M 2 ) 

—0 (by definition) 



(ill 



p 2 =M 2 



2G(p 2 - M 2 ) 



p 2 - M 2 



(Bl) 



(B2) 



, (B3) 



p 2 =M 2 



where in the last step we simply defined g\ Iq -q- By means of the wave function renormalization we define a normalized 
propagator in the following way 



s m (p 2 ) = gliagSuip 1 ) 



(B4) 



Now Sm(p 2 ) is a propagator for renormalized mass fields 4> = guqqfy- At this point we consider a complex mass pole 
solution 



P 2 = (M± 1 T -) 2 



(B5) 
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where for small V we get p 2 « M 2 ± iMr and the corresponding propagator, in a similar way, 

o ( 2 )= i = P 2 - M 2 ± zMr 

MlP J p 2 -M 2 TiMT {p 2 - M 2 ) 2 + (MT) 2 ' 1 j 

We introduce the spectral function by its definition to get 

~ Mr 

A(s = p 2 ) = 2* ImS M ~ ± (g _ M2)2 + (Mr)2 • (B7) 



